Effect of Disorder on Transport in Graphene 
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Quenched disorder in graphene is characterized by 5 constants and experiences the logarithmic 
renormalization even from the spatial scales smaller than the Fermi wavelength. We derive and solve 
renormalization group equations (RGEs) describing the system at such scales. At larger scales, we 
derive a non-linear supermatrix cr-model completely describing localization and crossovers between 
different ensembles. The parameters of this a -model are determined by the solutions of the RGEs. 
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Introduction- Dirac spectrum of the quasiparticles in 
graphene confirmed by recent experiments [1-3] is a con- 
sequence of the honeycomb lattice symmetry [4]. Al- 
though many properties of graphene can be understood 
in terms of ballistic motion of "relativistic" electrons de- 
scribed by a Dirac-like equation (see, e.g. [5, 6]), disorder 
plays an important role in sufficiently large samples. 

Influence of disorder on a two-dimensional electron 
gases on the honeycomb lattice were studied in several 
works [7-13] within a self-consistent Born approximation 
(SCBA) standard for weakly disordered metals and su- 
perconductors [14]. Within SCBA, such quantities as 
density of states, or localization-less conductivity [7-10] 
were calculated. The weak localization (WL) correction 
was discussed in Ref. [12] and calculated in Rcf. [13]. 

However, the SCBA is not justifiable for the Dirac 
spectrum, \e\ = v\p\, as it can be seen already in the 
fourth order perturbation theory in the disorder poten- 
tial, sec below, and a more careful analysis is needed. 

In this paper, we reveal the origin of the logarithmic 
effects specific for the Dirac spectrum and different from 
WL. These corrections, see Fig. 1, are contributed by all 
spatial scales between the lattice constant a, and either 
wave-length X e = hv/e, or the scale determined by disor- 
der, and that is why we will coin the name of "ultraviolet 
logarithmic corrections" (UvLC) for them. We will sum 
up the leading series of UvLC within a one-loop RG. 
At larger linear scales all the physics is described by a 
non-linear cr-model[15] and UvLC enter as renormalized 
parameters in this model. We will show that the low- 
energy asymptotics correspond to the orthogonal ensem- 
ble. Thus, the one-particle states are localized at any en- 
ergies in contradiction with the findings of Rcf. [1] of the 
minimal metallic conductivity in the undoped graphene. 

Disordered Hamiltonian- In the undoped graphene two 
bands cross the Fermi level at K and K' points. The 
corresponding Bloch functions comprise the basis of the 
four-dimensional (4d) representation of the planar group 
of the honeycomb lattice. We join them in a vector 



where we use the fact the points K and K' are connected 
to each other by the time reversal symmetry 



ip*(r) = Z(f(r); z 



KK' 



(2) 



where 4d-spacc of the wave-functions is represented as 
a direct product AB <E> KK', of the sublattice AB and 
"valley" KK' 2d-spaces, and t" z , t± = (t x ±it v )/2 are 
the Pauli matrices acting in those spaces, a = AB, KK' 
(we omit the physical spin). 





FIG. 1: Interfering scattering events involving (a) two- and 
(b) three-impurity scattering. They produce the correction 
(a) O[ln(X e /a )] and (b) O [in 2 (X e /a )]. 

The low energy properties are described in k ■ p- 
approximation [4], i.e. the wave function \E , e (r) is looked 
for as v P £ (r) = <^(r)-$ £ (r) with a smooth envelope, ^ e (r), 
satisfying the effective Schrodingcr equation (h = 1) 



Ho + V(r) $ e (r) = e$ £ (r), H = -ivEV, (3) 



where V = (d x , d y ) and we introduced 4x4 matrices [16] 
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fi T (r) = ((tpA, <Pb) ab ; (<Ps, -V*a)ab) 



AB) KK' 
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For the spinless particles, the time-reversal symmetry 
(TRS) requires (r) = <p* (r) • $* (r) also to be the eigen- 
state of the original Hamitltonian. Together with Eq. (2) 
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it constrains the effective Hamiltonian from Eq. (3) as 
H = zH^z, V = zV T z. (5) 



The disorder, V(r) = V^r), in Eq. (3) may break all 
the symmetries except of the TRS (5). This leads to 

V"(r) = lu (r) + Grm,iVm,i( T )> ( 6 ) 

where u , u rn i are real random functions, 1 = 1 AB ® 
\ KK , and the summation over indices {i, m} = {.t, y, z} 
is implied in Eqs. (6)-(7a). Equation (6) is the most gen- 
eral form for the disorder before averaging. Averaging 
must restore the rotation, reflection, (Cev), and trans- 
lation symmetries, which leads [17] to the most general 
Gaussian correlations (70 > 0, T m ^ > 0) 



(K ® v;) = S r 



(7a) 



Among the coefficients T m ^ only four are independent: 

r* = 7 *, rt*-»> = 7 ±, rk tf} = p„ r\ x x f } = p ± . (7b) 

In what follows, we will treat the five independent pa- 
rameters 7, j3 of Eqs. (7) as the starting point of the the- 
ory, not trying to calculate their particular values deter- 
mined by the details of the single impurity on the scale of 
the lattice constant. However, some qualitative conclu- 
sions can be drawn from the symmetry arguments. The 
diagonal term, (70 > 0,rj„ = 0), is generated by the di- 
lute impurities such that a single impurity in the crystal 
still preserves the point symmetry Cq v . Examples are 
the interstitial atom, or, more importantly, the remote 
charge impurity potential weakly varying on the scale of 
the lattice constant, a. On the other hand, a vacancy in 
otherwise perfect honeycomb lattice preserves only C^ v 
symmetry and allows for ~f z ~ 70. The intravalley scat- 
tering, 7x > arises e.g. due to the local deformation 
[18] or any bond disorder breaking the C3 symmetry. The 
latter also causes the intervalley scattering, (3 Z> ±. Finally, 
the terms of the similar structure are generated even by 
the multiple scattering by the pair of vacancies if the 
higher gradient terms are included in Hamiltonian (3). 

Origin of UvLC lies in the interference of the waves 
multiply scattered by the disorder potential, however, the 
configurations of the impurities giving rise to the effect 
differ from the familiar [15] WL. Thus, it is instructive to 
explain them, first, in terms of the counting of the scat- 
tering events and then present the rigorous calculation. 
It will also illuminate reasons for the failure of SCBA. 

If there were no disorder potential the wave function 
wou ld be the plane-wave with with a momentum k 
and the structure of Eq. (1). Limit ourselves by positive 
energies e = vk, k = |k| for concreteness, and the projec- 
tion operator P(n) = (1 + Sn) /2; n = k/|fc| chooses the 
chirality corresponding to the positive energy. Consider 



, u impurity placed at point Ri . The asymptotics of wave 
function acquire the form 4>k = $j t °' ) + ^ , wnere 
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and 17 = r — R4 is the distance from the impurity. 

Equation (8) is valid for both |ri| and X e much larger 
than the characteristic size of the scatterer. All the the 
details of the impurity are encoded in the matrix f\ which 
can be viewed as the scattering length, ao =|j fif\ H 1 ^ 2 - 
The value of ao is of the order of the size of the impurity 
which brings the estimate — a for neutral impurities. 

The asymptotic behavior of Eq. (8) at |ri| 3> X e is 
nothing but the out-going spherical wave, corresponding 
to the s-scattering. It is important to emphasize that 
the dominance of the s-channel is the consequence of the 
large wave-length and not of the peculiarities of the im- 
purity potential. The proportionality coefficient between 
the spherical wave and the amplitude of the plane-wave 
is the scattering amplitude that enables us to find the 
clastic scattering cross-section of the electron with the 
original momentum direction to the final direction xif 



ds 



(In 



f 



2tt\ e 



Pfr/J/iPfri) ® PCn^/iFfo) . (9) 



If there were no interference, Eq. (9) would describe 
all the kinetics of the system. For the impurity density, 
n-imp *C 1/og, the mean free path is estimated from 

iel - (n„ np II s Wy 1 ~ X £ / (n imp al) > X £i . (10) 

To understand the role of the multiple scattering, con- 
sider the two impurity scattering, see Fig. la(i)-(ii). 

Applying Eq. (8) twice we obtain for IR12I <C X e <C 
I r 1 1 , the outgoing wave with f\ — > 5f = /12 + /21 and 



Ry — R; — Rj, (11) 



for i, j = 1,2, j, and ao < IR12I < %e- Equation (11) 
describes the two-impurity scattering amplitude for the 
given configuration of the impurities. The transport, see 
Eqs. (9)- (10), is determined by the powers of the fuj\ 
averaged with respect to all configurations 
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where C = In (X e /a). Analogously, we find 

Mf{i3)®flji)) ~ ~ n «npC [fitjj] ® [fit a fi] . (12b) 

Combining Eqs. (12a)-(12b) we obtain: 
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f&afj ® flVJl , (12c) 
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where non- vanishing coefficients are d}?> 
d^f = — 1, and summation over all repeating indices here 
as well as over a = x, y in Eqs. (12), (14) is implied. 

Equation (12c) is the main result of the qualitative con- 
sideration revealing the origin of the logarithmic diver- 
gence. The Boltzmann kinetic equation systematically 
neglects those contributions. SCBA scheme accounts 
only for the diagonal components d\?, , d 2 \ and misses all 
the other contributions. For the scalar disorder, e.g. it 
leads to the violating the TRS of the problem. 

There are two more sources for the ni mp a 2 C correc- 
tions. One of them is the correlation of the one- impurity 
scattering with the two-impurity scattering in which one 
of the impurities is visited twice as shown on Fig. la(iii)- 
(iv). The corresponding result is easily obtained from the 
three- impurity scattering amplitude [cf, Eq. (11)] 

f m = (/,SR,,/,ER, ;../;) / (4 7 r 2 |R ij f|R ife | 2 ) , (13) 

and we find 



M/ (121) /t> 



'■imp 
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This correction is missing in SCBA. The last logarithmic 
effect arising in this order is the logarithmic dependence 
oc rii mp a 2 e In e of the averaged forward scattering ampli- 
tude. It does not affect scattering processes directly but 
renormalizes the spectrum of the free Hamiltonian Hq. 

Analogously, one considers the scattering from the 
three-impurity configurations, see Fig. lb), starting from 
the expression (13) and finds 54 contributions oc n 2 mp C 2 
only 6 of which are included in the SCBA. 

Interference of few processes, is somewhat special as 
they do not vanish completely even if the distance be- 
tween impurities is larger than X e . In fact the interfer- 
ence 1 — 2 and 3 — 5 are the first contributions giving rise 
to the WL correction which is also logarithmic. The WL, 
however, originates from the spatial scales larger than £ e i, 
and that is why they can be separated from the UvLC. 

Field theory and RG— Let us turn to the rigorous cal- 
culations using the supersymmetry method. Due to the 
4x4 matrix structure of the Hamiltonian, Eqs. (3, 6), 
the supervectors ip should have 4 times more compo- 
nents than usually used[15], i.e we need 32-components 
for calculation of the conductivity. The resulting 32d- 
space can be presented as a direct product of five 2d ones, 
AB (g> KK' ® AR <g> eh® g, where AR, eh and g are the 
retarded-advanced, particle-hole, and the fermion-boson 
sectors. Averaging over V and using Eq. (7a), we find 



■ cxp (-L [-0]) T>ip, ip^A = ip = 



Cip 



A = A <%> 1 AB : A = t ar ®1 kk ' 



i eh 



I 9 : 



(15) 
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where ... on the LHS stand for any combina- 
tion of advanced/retarded Green functions G A ' R = 



£Titi/2-JJ - VtM) and ••• in the RHS for 
the corresponding sources breaking the supersymmetry. 
Their form can be found in Ref. [15] but it will not be im- 
portant for us. The Lagrangian L [ip] = Lq [ip] + Li nt [ip] 
is given by (sum over i, m = x, y, z is implied) 



Lq [ip] = i J ip 

Lint VP] = o 



e - H - A ( - + iO 



tpdr, 



la {iPipf + T\ n (i>G m ^)' 



(16) 



dr 



i eh 



where (H , G m , 4 J = (H , G m , 4 J <g> I 9 <g> 1 AR <g> l e 

The perturbation theory in Li nt [ip] leads to UvLC and 
we can calculate the integral (15) using a RG scheme. 
We decompose ip as ip = ipo + ip, where ip is slow and 
ipo is fast, and integrate out ipo- We rescale ip and the 
coordinates to keep the coefficient in front of e and the 
ultraviolet cut-off intact. It gives back Eq. (16), with 
renormalized couplings 70, T\ n , Eq. (7b) and the velocity 
v, Eq. (3). This yields RGE (see Ref. [17] for details), 
which we display for the most interesting case 70 > T^: 



2-KvdtV = - (70 + g\\ + 2g±) ; 
9irv 2 d tl0 w 2 ( 5 f + 2ffi) ; 7rv 2 d t 5 9ll 
9irv 2 d m w -8gj - 20g l]g± + Ug 2 ± ; 
9-Kv 2 d t gj_ « 4g||gj_ - 18^. 



(17) 

-370^11 ,±; 



Here g\\ = -y z + 2-f±; g± = (3 Z + 2(3±; 6g\\ = j z - -f±: 
Sg±^ = (3 Z — /3± ; t is the logarithm of the running energy. 
Solving Eqs. (17) down to the energy, \e\, we find 



B (e)=|M 
7r e 



1/2 



7o (e) = 7o + C 



ff||( £ ) 
70 



(18) 



where £0 ~ J cxp (—^v 2 /jo) is the energy at which the 
first loop RG breaks down, and J is the bandwidth. 

The last of Eqs. (17) yields non-mixed valleys, g± = 0, 
to be unstable, and g± flows towards 

0||(e) « 9±{e) w 97o/{141n[r/ln| £ |/ £o ]} , (19) 

$9\\,±( £ ) m3 (| £ l/ £ o), and t* depends on g\\ t ±(e — J)- 

Equations (18) enables one to use the renormalized 
parameters for the standard calculation of the transport 
coefficients. The diffusion coefficient D(e) is given by 



D(e) = v 2 (e)T tr {e)/2; l/r tr (e) = ir l0 is e /4, 



(20) 



where v e = \e \ / (ttv 2 ) is the density of states (per one 
physical spin). Einstein relation and Eqs. (18) yield [20] 

4e 2 



n 2 h 



ln(H 



(21) 



Equation (21) is the universal formula for the UvLC. 
It describes either the temperature or the density de- 
pendence of the conductivity limited by the short-range 
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disordered potential |e| — > max (ef,T). It also gives 
the leading dependence of the thermopower through the 
Mott-Cutler formula. At |e| < e , the logarithms are cut 
by 1/Tt r , which leads to the replacement of ln(-) to the 
factor of the order of unity. The precise value of this 
factor, however, can not be calculated within Eqs. (17). 

Nla-model and localization- At distances larger than 
vj max(|e|, £o) soft modes giving rise to UvLC freeze out, 
and only the degrees of freedom guarded by the funda- 
mental symmetries of the system remain gapless. All 
those degrees of freedom are described by nicr-model [15] 
and functional integral (15) is replaced by integral over 
16 x 16 supermatrices in KK' (g> AR(g> eh® g space [17]: 

(...)= J ■ ■ ■ exp (-F [q] ) T>Q, Q 2 =t- Q = CQ T C T ; 

(22) 



0(T) 
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The free energy functional F [Q] takes the form 



2iujAQ 



y z9\\ 



Px,Q 



Pv,Q 




(p a ,%) = ( 



t5 K '<8> t eh ,t KK '®T eh 



) ® t AR ®t 9 , (23) 



where A is the vector potential due to the magnetic field 
normal to the graphene, and the entries include UvLC. 

Equation (23) is the only form allowed by the symme- 
tries of the problem. The symmetries of the Q-matrices, 
Eq. (22), correspond to the two replicas of symplcctic 
ensemble, which would flow to the limit of large con- 
ductances. However, due to gu j_ > only Q cx t KK 
is allowed at large distances and one obtains a generic 
orthogonal ensemble. Thus, all the eigenstates are lo- 
calized. Schematic temperature dependence of a for un- 
doped graphene is sketched on Fig. 2, [22]. 

The fist loop correction in Eq. (23) yields the WL [20] 
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where -d = di >2 ,3 = 1, \]_ = ZxVe9\\,±, 1 b 
4D(e)eB/(hc), Y(x) = ip(x + 1/2) +lnx and ip(x) is the 
di-gamma function. Inelastic processes are accounted for 
by — iuj — > — iuj + t^ 1 , where is the dephasing time. 

Equation (24) agrees with Ref. [13]. The new informa- 
tion here is the logarithmic dependence of the parameters 
on the electron energy [21], see Eqs. (18)-(20). 

In conclusion, we have presented a complete descrip- 
tion of a disordered graphene and demonstrated that 



(24) 



Weak localization, Eq. (24) 
Orthogonal Symplectic 



UvLC, 
Eq.(21) 



T~E n 



In T 



FIG. 2: Schematic dependence, <r(T), for the undoped 
graphene and for r^ 1 ocT. 



there are two different types of logarithmic contributions 
into physical quantities, see Eqs. (21) and (24). 
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